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Os ! Abstract 
On 

£>V We show that all abnormal parity three-body meson interactions can be ade- 

quately described by quark loops, evaluated at zero external momentum, with 
couplings determined by U(Nf) symmetry. We focus primarily on radiative 
meson decays which involve one pseudoscalar. The agreement with experi- 
ment for non-rare decays is surprisingly good and requires very few param- 
eters, namely the coupling constants g wqq and g pqq and some mixing angles. 
This agreement extends to some three-body decays that are dominated by 
■ pion pairs in a P-wave state. 

in 
in 
o 

G\ 

^ ■ One can perfectly well view tt° — > 77 decay as arising from the PVV quark-level JIJ 

linear a model (LcrM) instead of through the AVV chiral anomaly . In such a description 
Ph. it is straightforward to predict many other PVV radiative decays through quark triangle 

graphs. The purpose of this paper is to investigate how well this quark loop picture accords 



X 



I. INTRODUCTION 



with experimental data [[|; in this way we avoid discussing anomaly predictions of these 
processes and notions of strong CP violation with consequential (low mass) axions, since 
these effects have yet to be discovered in the laboratory. We also consider the situation for 
rare PVV decays that are associated with tiny quark mixings (such as J /if) — *> np, r)'<p) 
and with isospin-changing corrections (such as J /if) — > ttuj, rjp). 

In Section II we show that the dominant LaM quark triangle PVV predictions for 
P — > 77, P — > I/7 and V — > -P7 all fit the data, within experimental errors; the only 
parameters at our disposal are gp qq and gv qq , the couplings of the pseudoscalar and vector 
mesons respectively to the quarks, the constituent quark masses m q and a couple of well- 
known mixing angles. We pursue this approach in Section III, for the rarer isospin-conserving 
interactions where the data gives valuable information about the small amount of admixture 
of light quarks in the heavy qq states, and vice versa; here the agreement is less good, but 
the experimental information is also sparser and may be expected to change with time. In 
section IV, we show that a simple picture of AI = 1 meson-meson transitions, associated 
with the u — d quark mass difference, nicely provides the rates for isospin changing decays. 
The last section contains some preliminary work on three-body decays in which at last one 
pair of pions is produced in a P-state; we show that p-meson poles essentially saturate 
the amplitudes for those processes, given the two-body amplitudes determined previously. 
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The reasonable agreement between naive theory and experiment suggests that this quark 
loop picture is a viable alternative to the field-theoretic anomaly-inspired evaluation of such 
amplitudes. However the tantalizing picture which emerges is that internal mass dependence 
of the box diagrams is largely negated by the form factors which must be present at the 
meson-quark vertices. Such are our conclusions in Section VI. 



II. QUARK TRIANGLE GRAPHS FOR PVV INTERACTIONS 

For the past thirty years, the tt° — > 77 decay has been understood as an AVV chiral 
anomaly 0] combined with the partial conservation of axial current (PCAC). Writing this 
process and all other PVV amplitudes in the conventional form Mpyye^^^k' 1 'e K k* ', where 
k and e refer to vector particle momentum and polarization respectively, one derives 

M n0jJ = -f^L- = ^ « 0.025 GeV" 1 . (1) 

With f n pa 93 MeV and a quark color number N c = 3, the theoretical result is accurately 
consistent with data H: 



I M^ 7 |= y , 64vrr w o 77 /m3 = 0.025 ± 0.001 GeV . (2) 

The question then arises as to which theoretical models also predict the same result while 
simultaneously demanding that N c = 3. 

One such favored model is the quark-level jl| linear a model (LerM). This model also pro- 
vides u and d quark triangle graphs but with the meson interacting through the pseudoscalar 
7 5 vertex, rather than through the divergence of an axial vector. No accompanying LaM 
it + loop is allowed because the 7r°7r + 7r~ coupling vanishes; consequently the AVV anomaly 
result immediately follows in the spirit of Steinberger 0j . (In carrying out the evaluation of 
M, we go to the soft limit, because the pion is essentially massless, relative to other hadronic 
masses and even the constituent quarks.) The fact that this LaM also demands N c — 3 |1J 
is tied to the Lee null tadpole condition || and ensures dynamical generation of masses, in 
contrast to the original nucleon- level LaM |J. The pion-quark coupling g nqq needed obeys 
the quark-level Goldberger-Treiman relation (GTR) g nqq = m ns / f n where the constituent 
(nonstrange) quark mass m ns = rh is dynamically generated to the expected M N /3 value in 
the chiral limit, since N c = 3 is required now. 

Our aim in this paper is to investigate how well the sigma model picture fits the data for 
PVV processes. We take the naive viewpoint that the PVV amplitudes, just like 7r°77 can 
be estimated by working out Mpw i n the limit of zero external momentum, even though 
this is a long way from the mass shell for the heavy mesons particularly. The reason why 
we do so is because many of these processes are quite close to particle thresholds and would 
otherwise be highly dependent on binding energies and the like — for which there is no 
evidence. In such a soft limit and including the crossed diagram, it is very simple to work 
out that 

w 9pi9v 2 9v 3 N c rm 1 a + m 2 l3 + m^ 

Mp ^ = 2** J mla + m^ + mh 6{1 " " " fi " l] ^ ^ (3) 
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where rrii are the internal quark masses, in an obvious cyclic notation. The Feynman para- 
metric integral is especially easy to work out when all quark masses are equal to m, when 

w 9Pi9v 2 9v3 N c 

It is not much harder to determine the integral (3) if two masses are equal, say mi = m 2 = 
m, m 3 = m', when 

9pi9v29v 3 N c J{r) _ 2r 4rlnr _ m 

2n 2 m + m r 2 — 1 (r 2 — l) 2 m 

Thus J(l) = 2 in the equal mass limit m = m'. See Figure 1 to see how J(r) varies with r. 

It only remains to substitute the couplings of the mesons to the quarks in order to derive 
the theoretical values for M PV v- To do so we will adopt the naive perspective that all strong 
couplings are related by U(Nf) symmetry. Therefore, in the absence of any isospin breaking 
and neglecting the third generation of quarks, let us see how far one can go with the effective 
Lagrangian (suppressing the 75 matrix between spinors), 

£pqq — gp[^°(uu — dd)j\pl + n + ud + Tr~du + rj ns {uu + dd) / y/2 
+K + us + K'su + K°ds + K°sd + r] s ss + r] c cc 
+D°uc + D~dc + D°cu + D + cd + D;sc + Djcs]. (4) 

The universal coupling gp above is fixed by the 7r° coupling, which in turn is determined 
by GTR; thus g P = (m u + m d )/y/2f^ = \f2ihj f w ~ 5.13. We must also relate the physical 
neutral pseudoscalar states to the isospin combinations in Eq. (4) via certain mixings. 
For the purpose of this and the next section, we shall disregard electromagnetic effects in 
the strong sector and take (see Appendix) 

Vns = V cos (f) P + 7]' sin (f) P ; r] s = -r]sm(f)p + r/cos0p; (f) P ~ 42°, 

with i] c unmixed. (Later on we will examine rare decays, where the small residual mixings 
become rather important.) The strong vector meson interactions are likewise written as 

£>Vqq = 9v[p°(uu — dd)/V2 + p + ud + p~du + u ns {uu + dd)/V2 
+K* + us + K*~su + K*°ds + K*°sd + uj s ss + u c cc 
+D*°uc + D*~dc + D*°cu + D* + cd + D*-sc + D* s + cs], (5) 

suppressing vector indices and the 7 matrix between quark spinors. Here the coupling 
constant is fixed by that of the p-meson which is itself determined by the leptonic decay 
rate: g v = g p /\/2 ~ 3.56. These vector states also undergo the semistrong mixing via 

u) ns = ucosifiv + (psin(p v ; u s = — u sin0y + 0cos</v; 0y — 3.8°. 

(At this level we may identify u c with the J/ip particle.) Finally, it is trivial to consider 
radiative decays, by substituting 7 for the meson and including the proper electromagnetic 
coupling of the photon to the quarks, e q . 
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The magnitudes of the Mpyy can be found directly from experimental decay rates, as 
they occur ||. For P — > VV decays, the rate is given by 

Tp VV = A 3 \M PVV \ 2 /32iTm P = p 3 vv \M PVV | 2 /4tt, (6) 

where A = 2mppyy and pyy is the magnitude of the three-momentum of one of the vector 
mesons in the rest frame of the decaying P. (We must be careful to multiply the right-hand- 
side of (6) by a factor of 1/2 when the two vector mesons are identical.) A very similar 
formula applies to the decay rate V — > PV, 

T VPV = A 3 \M PVV \ 2 /967rm 3 / = p 3 VP \M PVV | 2 /12tt, (7) 

where my refers to the mass of the initial vector and pvp = A/2my is the momentum of one 
final particle in the decay rest frame. It is then a simple matter to extract the \M\ from the 
measured rates T, for subsequent comparison with theory. We shall do this constantly, 
without further elaboration. 

We mentioned at the start that the M n o 77 amplitude is theoretically determined in the 
sigma model by ing the u and d loop contributions and equals e 2 gp/Ay/2n 2 m = a/nf^, 
fitting experiment admirably. Indeed many of the processes involve neutral mesons and 
require the interaction terms, 

gp[ir°(uu — dd)/V2 + (rj cos^p + rf sin (f)p){uu + dd)/V2 + (—rj sin <\> P + rf cos 4>p)ss + r] c cc], 

(8) 

gy[p°(uu — dd)/v2 + [u cos<py + (p sin <py){uu + dd)/V2 + (—uj sin0y + <p cos 4>y)ss + i/jcc]. 

(9) 

Instead of treating each process one at a time, let us consider just three examples to explain 
what we are doing, before presenting our results for non-rare decays in a table. 

Consider first the decay rj' —* w-f, which contains a charge coupling, one vector coupling 
and one scalar coupling. The sum of the loops of nonstrange quarks and the mixing angles 
implied by (8) and (9) gives (e = V&m ~ 0.3028) 

Mrf^ = eN c g v g P sin (f) P cos 0y/247r 2 m ~ V^egy sin (j)p/8iT 2 f^. (10) 

Using the values quoted previously, we obtain the rough theoretical magnitude, M n > wi ~ 
0.139 GeV -1 , which can be compared with the experimental result, 

= yJ^T^/pty * 0.137 GeV- 1 . 

Our second example is the semistrong process — > 7rp, which is governed by nonstrange 
quark loops but is rather sensitive to the vector meson mixing angle <py. The theoretical 
result here is 

M^o = V2N c g v g P sin <fi v /8ir 2 m ~ 3gy sin ^y/Air 2 ^. (11) 

Inserting the accepted value of (fiy ~ 3.8°, the prediction is that M^o p ~ 0.69 GeV^ 1 , and 
this is in agreement with the upper bound obtained by experiment: M^ p < 1.18 GeV -1 ; 
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in this connection we note that the 1996 PDG compilation cites a specific value for 
this particular channel (which we prefer), but this has been retracted in the 1998 PDG 
compilation, resulting in less predictive power. Finally we take a look at the decay ip — > ry c 7, 
which is governed by a charm loop and involves one photon. This time we get 

Mihci = eN c g v g P /6TT 2 m c ~ eg v m/ V2tt 2 ferric, (12) 

which we estimate to be about 0.19 GeV -1 , for a constituent quark mass ratio of m/m c ~ 
338/1500 ~ 0.225. This is close to average experimental value of M^l p ~ 0.17 GeV -1 , 
obtained from the decay rate. 

In this manner we can work out theoretical estimates of all dominant PVV processes. 
We have tabulated the results and the experimental magnitudes for comparison in Table 
1. The near agreement of naive theory with data is actually quite astonishing, because we 
are dealing with sizeable external masses in many instances, particularly for heavy charm 
quarks, and so are a very long way from the chiral limit. We will return to this point towards 
the end of the paper. 

It should be noted that the last five decay channels in Table I correspond to unequal 
mass quarks running round the loop, in contrast to all the previous cases, and the absolute 
value of the amplitudes are perhaps not as impressively predicted are the others (although 
the ratio \Mx*+K+-yl / M^»o^o 7 | is reasonably close to experiment). On the whole the fit is 
quite satisfactory, considering the fact that we have no free parameters at our disposal. 

III. RARE ISOSPIN-CONSERVING DECAYS 

There are many processes involving the heavier quarks which conserve isospin. They 
tend to be rarer than the cases considered above and measure the smaller admixture of 
lighter quarks in the meson composites. In fact these rare processes provide much valuable 
information about such small mixings and, emboldened by the success of the soft momentum 
predictions of PVV amplitudes, we shall apply the same method to probe the mixing angles 
of neutral vectors and pseudoscalars. Let us write the meson which couples to the cc as the 
combination r] c + 5 c r] + 8' c rj', where the two 8 are very small quantities (which is why we have 
not modified the normalization factor in front of r] c in first approximation). It follows then 
that the non-strange and strange pseudoscalars are (in the same limit), 

Vns = rj cos (pp + rf sin (f)p — r] c (5 c cos (ftp + 8' c sm(j)p) 
t] s = — rj sin <p P + rj' cos <f) P + r) c (5 c sm <f)p — 5' c cos<f)p). (13) 

Consequently the interactions (8) of the neutral pseudoscalars get modified to 

gp[7r°(uu — dd)/\/2 + r]{cos(j)p(uu + dd) /y/2 — sin (fipss + 8 c cc} 
+ r/{sin cpp(uu + dd)/y2 + cos (ppss + 8' c cc} 
+ r] c { — (5 c cos (f)p + 5' c sin <f>p)(uu + dd)/\/2 + (<5 c sin <\>p — S' c cos 4>p)ss + cc}]. (14) 

Much the same idea can be applied to the neutral vector interactions (with mixing param- 
eters 8 replaced by e): 
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gv[p°(uu — dd)/ ' y/2 + ui{cos(j>y{uu + dd)/y/2 — sin (f) V ss + e c cc} 
+ 0{sin 4>y(uu + dd) / V2 + cos (pyss + e' c cc] 
+ ip{ — (e c cos0y + e' c sin0y)(-uM + dd)j \[2 + (e c sin 0y _ e ' c cos0v)ss + cc}]. (15) 

We are now in a position to tackle rare PVV processes involving r\ c and ip. These are 
evaluated in the same way as before, but now using the neutral couplings (14) and (15). The 
results are summarized in Table II; theoretical values are only worked out to first order in 
5 or e, since these are already small quantities. We need to give some words of explanation 
as to what entry in the last column is being predicted and what is not, in contrast to Table 
I, which is essentially parameter-free. In our parametrizations (14) and (15), the cleanest 
predictor of the amount of nonstrange and strange quark components in r\ c are the two 
decays r\ c — > 2p and rj c — > 20; they determine separate combinations of the mixings 8 C and 
8' c (which represent the amount of charm quarks contained in 77 and rj'). We may derive 
a reasonable fit by choosing the values 5 C ~ 0.0054 and 5' c ~ —0.0008; this permits us to 
predict the decay rates for r\ c — > 2u and r\ c — > K*K* in the table. While the latter prediction 
is good, the former is decidedly not - in fact the experimental value is slightly embarrassing 
and difficult to understand on any sensible theoretical basis; for U(2) symmetry between u 
and d entails that the process i] c p°p° should be almost equal to i] c ujuj, which is a far cry from 
what seems to be observed! 

With respect to the e parameters (because (f>v is so small) the process ip — > 7r°p° tells us 
about the nonstrange-charm mixing e c almost at once, while ip — > r}4> provides the mixing 
between strange and charm states, namely e' c . We find a fair fit (including sign) with the 
values, e c ~ —0.0002, e' c ~ 0.00009. This enables us to predict the amplitudes for the 
other entries in Table II involving e. They are all roughly correct, except for ip — + 7/7 
and ip v' 00 - We do not understand the reasons for this; blaming the discrepancy on the 
effect of mass extrapolation from the shell to zero-momentum would undermine the other 
reasonable answers; but in any case it should be noted that the vector meson dominance 
(VMD) prediction is a factor of 250 larger than the data which heightens our suspicions. It 
is possible that the experimental results may shift a little over time and bring the results 
into closer agreement with theoretical expectations, or possibly the experimental decay rate 
also includes the contribution from a (980). 



IV. RARE AI = 1 VPV DECAYS 

Here we want to take a look at four processes, 

1p — > 7T°CJ, lfj — > 7T°0, 1p — > 1]p°, 1p — > Tj'p , 

because they involve isospin breaking interactions of mesons, without being accompanied by 
photons. (Interestingly, there seems to be no comparable data in the r] c decay sector.) All 
four processes are smaller by at least a factor of three compared to processes that conserve 

isospin; eg MJX - - 3M Sv>- 

We may ascribe the existence of these small amplitudes to electromagnetic effects and 
the u — d mass difference. With regard to the quark triangle loop, the u — d difference gives 
a correction of order ~ 5/350, which is insufficient to explain the observed magnitudes. We 
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shall find below that the experimental results can be reasonably well obtained by AJ = 1 
mixings between mesons (including the well-known p — uo transition), which are themselves 
dominated by the u — d difference in quark-loop self-energy contributions. 

Before going further, let us estimate two strong PVV amplitudes that cannot be directly 
measured from decays, because of the masses of the participating mesons; they are nonethe- 
less important in a subsidiary role for what follows. Firstly, there is the interaction p°ir°uj 
(which actually enters uj — > 3tt, via pole dominance), that we estimate to be 

Mpo^ = zV2glg P cos 2 (p v /8TT 2 m ~ 3g 2 p /87r 2 f n = 3C ~ 10.4 GeV" 1 . (16) 

Next we shall need the amplitude 

M p o vp0 = 3Ccos0 P ~ 7.85 GeV _1 . (17) 

Finally we shall require amplitudes which we have estimated previously fairly well (or which 
we can take directly from experiment), namely M^o p o ~ 0.0021 GeV -1 , M^ vw ~ 0.0016 
GeV" 1 and M^o p o ~ 0.69 GeV" 1 . 

Other quantities we will need are the (P\H\P') and (Vl-H^V) elements, where unprimed 
states refer to / = composites while primed states correspond to / = 1 mesons. These 
elements include the most well-studied case of the p — u transition, which is the sum of two 
pieces ||: 

= (p\H AI=1 \u) = (p\H JJ \u) + (p\H u - d \u) 

~ 700 - 5100 ~ -4400 MeV 2 (18) 

In similar vein, we may estimate the p — <fi transition to be 

M p4> = (p\H\(p) ~ tan(p v (p\H\u) ~ -290 MeV 2 . (19) 

We must be a little more careful with the p — ip transition, which is dominated by the photon 
exchange contribution H JJ \ 

M P1P = <p|tf A/ =V> = (p\H JJ \^) + (p\H u - d \^) 

~ (em p m^) 2 /g p g^ + e c (p\H u ~ d \u) 

~ 1010 + 1 ~ 1010 MeV 2 . (20) 

Last, but not least, we shall require the n — r] and 7r° — transitions. This arises from a 
quark self-energy bubble plus an a tadpole term. The latter equals the former contribution 
U because g aoms n = (^ 2 - ™% na )/f n by SU(6) and chiral symmetry (JDJ. Noting the gap 
equation Q stemming from the GTR, the total result sums to 
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(27r) 4 p 2 — m 2 d p 2 — 



m 



2 - 



~ 2(m 2 - ml) ~ -5400 MeV 2 for m d - m u ~ 4 MeV. (21) 
We thereby estimate the two neutral pseudoscalar elements, 
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M m = (r]\H\n°) ~ -5400 cos P ~ -4000 MeV 2 . 

M v > n = (r]'\H\n°) ~ -5400 sin P ~ -3600 MeV 2 . (22) 

Actually the last element leads to very small corrections because each of the mesons is rather 
far from the mass shell of the other. 

We are now in a position to estimate these rare amplitudes. (Remember that we may 
ignore the negligible correction induced by u — d mass differences within the triangle loop.) 
The first case is 

M^ u = M t Mp °? + M T M "T + - 0.0007 GeV -1 , (23) 

m^p — nip — m l p — mi 

upon substituting the AI = 1 transition elements found above and the vertex amplitudes 
determined previously. Similarly, we find 

= M *' M W + »f + M^M±± ~ 0.0009 GeV 1 , (24) 
m^ — nip — m l p — 

M*i,„M n o„ n o M n ,,Mj mil , M mr M,i„ojo „ i , 

AW = f pvp + T T + v Z - 9 P ^ 0.0005 GeV" 1 , 25 
m^, — m~ — m~j m~ — m% 

M Wp o = % p pnp + p " ^ + ^ ~ 0.0006 GeV" 1 . 26 

— Trip nip — — m* 

The experimental values of these amplitudes, deduced from the measured decay rates, are, 
in units of GeV -1 , 

I m ;XI = 0.00067 ± 0.00004, |M^| < 0.0009, 

\M^fl\ = 0.00048 ± 0.00003, \M™f p0 \ = 0.00040 ± 0.00004. (27) 

All of these answers are in the "right ball park" and we venture to conclude that this A J = 1 , 
pole-dominated mechanism captures the main features of such feeble decays. 



V. SOME THREE-BODY DECAYS 

Here we shall have a look at a few three-body decays, in which pions are produced in 
a P-state and are known to be dominated by p-meson poles; specifically we shall examine 
the three processes ui — > 3tt and rj, rf — > 2tt^ to see how well they tie up with the earlier M 
amplitudes. Before doing so we shall require more accurate values of the on-shell couplings 
g P wiT, 9wp-K- We recall that the quark loop alone gives g pn7T = V2gv — 5.03 and g W(m = 
3gy I '47T 2 f n ~ 10.3 MeV -1 . However, in a sigma model, besides the non-strange quark loop 
we have to add a meson loop associated with sigma exchange. The effect of this is to enhance 
the value of g p7r7r by a factor of 6/5, to about 6.1, which is very close to the experimental 
value coming from the p decay width. Correspondingly, the p-oj-rt coupling is enhanced 
§|TT[ by a factor of (6/5) 2 to 3g%j87T 2 U ^ 15 GeV -1 . 
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The amplitude for the VPPP process lu — > 3tt may be written in the covariant form 
^4w7ri7r 2 7r3 = ^KXp.v^VxV^V^^-bi-K-K-Ky where the scalar amplitude M is dominated by p mesons 
T3] in each of the three possible two-body channels: 



+ : + 



m 2 p — s m 2 — t m 2 — u 



1480 GeV , (28) 



In deriving this result, we have averaged over the Mandelstam variables: (s) = (t) = (u) 
m t/3 + — 0.223 GeV 2 . If we then follow Thews' phase space analysis [13], we get 



T(u -> 3tt) = \M m \ 2 mlYj(2iif (29) 

with the constant matrix element giving Y w = 4.57 x 10~ 6 . In this way, we predict T(uj — > 
3/t) ~ 7.3 MeV, very close to the observed rate of 7.5 ± 0.1 MeV. 

Our next object of study, rj — ► 7r + 7r~7, is dominated by a p° meson pole, but just in the 



two-pion channel [[L4j]. The process also takes the same covariant form as the previous case. 
Here we estimate 

M w ^ = 2g p7VW M rm /(m 2 p - s) ~ 9.74 GeV~ 3 , (30) 

upon substituting the rj-p-j amplitude found earlier (Table I) and going to the soft pion 
limit. Once more, we find 

Y( v -> 7T7T7) = \M riW7rl \ 2 mlY v /(2n) 3 ~ 55 eV, (31) 

for a phase space factor of Y v = 0.98 x 10~ 5 . The prediction (31) compares well with 
the observed rate of 56 ± 5 eV. 

Finally we consider the r/ — > 7r + 7r _ 7 decay. Unlike the previous two examples, the p° pole 
is now well within the physical region. The latest particle data tables state that T(i]' — > 27r7, 
including nonresonant contributions) = 61 ± 5 keV, whereas our prediction (see Table I) is 
that T(t/ — > pj) = |M )7 /p 7 | 2 p 3 7 /47r ~ 65 keV. This strongly suggests that any nonresonant 
P-wave contact contributions are quite small and that the p poles accurately match the data. 
Incidentally it also confirms that the pseudoscalar mixing angle is <pp ~ 42°, rather than 
<p ~ 35°, coupled with the anomaly, as the latter angle gives a rate for rf — > titt^ which is a 



factor of 20 less than the observed rate 15 



VI. CONCLUSIONS 

Taking stock of the results, we appear to have succeeded in estimating some 34 VVP 
amplitudes with reasonable accuracy (see the Tables) barring two or three rare decays in- 
volving 7/ and uj mesons — which are themselves extracted from one or two experiments that 
surely deserve repeating in the future. The success extended to a few abnormal parity three 
body decays which are known to be dominated by P-wave pion pairs. We do not regard the 
coupling constants gp and gy which provide the overall scales to be 'parameters', because 
they are predicted in [1] or are derived from other processes and have well-established mag- 
nitudes; nor do we think that the quark masses are variables at our disposal, since they are 
also fixed by other means. The mixing angles <pp and (f)y are parameters but they too can be 
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found through processes other than the ones which we have considered in this paper. The 
only (four) variables which we were truly able to adjust were the mixing angles describing 
the non-cc content of i] c and ip. 

Perhaps the most surprising aspect of the analysis is the fact that the amplitudes, which 
we estimated in the soft limit, have reasonable magnitudes even for external heavy meson 
cc composites. We really did not expect this and we think that the conclusions are rather 
deep and point to some underlying heavy mass scale, way beyond the usual QCD scale of 
about 300 MeV. To put this deduction into perspective, consider some field-theoretic model 
which includes meson interactions with spinors; assume it is renormalizable for definiteness. 
Then the equation for the proper vertex function T(p, q) (where p is the meson momentum 
and q is the relative fermion momentum) will take the form 

r(p, q) = Z 1 + g 2 J d 4 q' S(p/2 + q')T(p, q')S(-p/2 + q') * K p (q, q'), 

where K is the appropriate interaction kernel for the model in question, S stands for the 
spinor propagator and Z7 is the renormalized contact term. If one specialises to pseudoscalar 
mesons and ignores fermion dressing for simplicity, the equation has the generic scalar form, 

A(p, q) = Z p + g 2 J d A q' A(p, q') * K p (q, q')/[p 2 /4 - q' 2 - m 2 }. 

The subscript p has been attached to the renormalization constant Z to remind ourselves 
that the renormalization procedure is often momentum dependent, although the infinite part 
of Z (or a l/(D — 4) pole term in the dimensional method) is of course insensitive to p. The 
limit as p — > of the equation above is readily taken and describes the meson vertex in the 
(off-shell) soft limit. 

If the pseudoscalar field is actually a composite of the fermions, then Z must vanish on 
the meson mass shell (p 2 = M 2 ) in order that the vertex equation reduces to a homogeneous 
Bethe-Salpeter equation. The results which we found indicate that the products of vertex 
functions and quark propagators running round the loop are rather insensitive to external 
momentum since they suggest that an amplitude like A(p, q) above depends very little on 
p, i.e. the g-dependence on the right hand side has compensating effects from form factors 
and propagators. (It is striking that the propagator denominators would appear to depend 
sensitively on the binding energy near the constituent threshold p 2 = 4m 2 in the infrared 
region of q', and there is no sign of this.) Therefore it indicates that Z p depends very little 
on p and that the homogeneous Bethe-Salpeter equation can very likely be extrapolated 
to p = with relatively little change. Our guess is that such loop integrals are really 
dominated by the ultraviolet region and a mass scale A of at least 10 GeV — much higher 
than the standard QCD scale — associated with the ratio p 2 /A 2 . Indeed, the seemingly weak 
dependence of our results on p are the most extraordinary part of this work and trying to 
understand the issue properly is an exciting avenue of future research; our discussion above 
does scant justice to the problem. 

ACKNOWLEDGMENTS 

This research was supported by the Australian Research Council. MDS appreciates the 
hospitality of the University of Tasmania where a substantial proportion of the research 
work was carried out. 



10 



APPENDIX: PSEUDOSCALAR MIXING ANGLE 



Since we are using quark loops to predict PVV amplitudes involving 77 or rf pseudoscalar 
mesons, it is imperative that we determine the most important 77 — 7/ mixing angle <pp from 



processes other [16j than PVV decays. Here we shall consider two other ways of fixing this 
angle. 

First we consider tensor (T) — > PP decays, a 2 (1320) — > KK, rjir, 7/71". Given the 
respective percentage branching ratios, from the PDG tables, 4.9 ±0.8, 14.5 ±1.2, 0.53 ±0.09, 
and remembering the p 5 phase space, we may deduce the two independent ratios, 



B(a 2 -> rpr/KK) = (prm/Px) 5 ■ 2 cos 



2.96 ±0.53 



(Al) 



B(a 2 -> t/tt/V) = (p^/iv) 5 tan2 0p = °- 037 ± °- 007 - 



(A2) 



Inserting the known magnitudes {p^iP^'iriPK) = (535,287,437) MeV, into this pair of 
equations, we may conclude from the first that <pp = 43 ± 5°, and from the second that 

Ki]/Ktt branching ratio of 0.003 requires [17 



42 ± 3°. Likewise the measured K% 



(pP = 41 ± 4°. As well, there are the tensor decays /, /' — > 7r7r, KK, 7777, but then the 
tensor angle <p T enters the analysis; yet <pp ~ 42° survives the results of the analysis. 

A second way of extracting 0p is through the QCD quark-annihilation diagram, involving 
at least two gluon exchange. Diagonalization of the NS-S meson mass matrix leads [IB] to 
the f] — rj mixing angle, 



arctan 



(m 2 , — 2m 2 K + ml)(m 



{2m 2 K — m 2 — m 2 )(m 2 , 



mi 



1/2 



41.9°. 



(A3) 



The close agreement between <frp determined by T — > PP data and QCD theory justifies 
our use of <pp ~ 42° in Tables I and II. 

In the singlet-octet basis, the mixing angle 9 is instead given by fl7| , |T8[ 



6 = 6 — arctan \pl 



54. T 



or 



cos 1 



—= cos 6 
\/3 




2 1 

sin <t> = \ — cos 6 H = sin 6. 

^ V3 



Thus for the pseudoscalar mesons, 6p ~ —13°, which is midway between the original Gell- 
Mann-Okubo value of —10° and most recent determinations of —20°. In this connection, it 
is interesting to compare the quark triangle predictions, 



\M, 



a 



2m 



(5 - — ) cos 9 P - \/2(5 ± — ) sin 9 

m R 



(A4a) 



I ~^^~t^ t yy I 



a 



(5 ) sin 6 P + V2(5 H ) cosflp 



with the predictions obtained via the AVV anomaly [H| . 



(A4b) 
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I i\janom I 



—j= ^j- cosOp — 2v / 2 ; y- sin^p 

V37lf n [fs JO 



(A5a) 



I t\ it anom I 
1^77 I 



V3nf n [fs Jo 



Comparing (A. 4) with (A. 5), we see that 

h/U = 3/(5 - 2m/m s ) 



0.83, 



(A5b) 



(A6a) 



fo/U = 6/(5 + m/m s ) ~ 1.05. (A6b) 

Both (A.6ab) exhibit U(3) symmetry in the limit m s — > m. However for the actual con- 
stituent quark mass ratio (obtained from the GTR or magnetic moments) of m s /m = 
{2Jk/ f-n — 1) — 1.44, the numerical values in (A. 6) follow. Note that the anomaly ver- 
sion |15j estimates fs/ — 1-3 via chiral perturbation theory but also finds fo/ f n — 1.04, 
which is close to our prediction; we agree with their justification that f / f n being near unity 
is what one expects if the singlet state and the pion have the same wave function. 
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TABLES 



TABLE I. Dominant strong and radiative decays that are somewhat insensitive to mix- 
ing angles. Under the column of theoretical formulae, Mpyy, we use the abbreviations, 
A = e 2 / Air 2 f 7T = a/irf 7T ~ 0.025, B = eg v g P /8ir 2 m ~ 0.207, C = g 2 , / Air 2 f n = V2g 2 v g P /S^ 2 m ~ 3.47, 
in units of GeV -1 , and r q = rh/m q for the constituent quark mass ratio. The latter are determined 
by the assumed values (in GeV) rh = 0.34, m s = 0.49, m c = 1.5, mj, = 4.8 



Process 


1 c ^(MeV) 


\M v \ (GeV J 


Mpw 


M y (GeV 


o 

7r — > 77 


/" t ^ /I 1 n ro\ ^, in — 6 

(7.74 ± 0.52) x 10 


n no r i i n nno 

0.0251 ± 0.008 


a 

A 


n no c n 

0.0250 


t? -> 77 


(4.6 ± 0.4) X 10 


0.024 ± 0.001 


/i[5cos (f)p — y2r s sm<pp\/6 


0.0255 


r? -> 77 


(a o i n a\ 1 n — 

(4.3 ± 0.4) x 10 


n noo i n nm 

0.032 ± 0.001 


Air' l i /o i 1 / o 

^4 [5 sin <pp + y2r s cos <ppJ/3 


0.0335 


% 77 


(7.5 ± 2.5) x 10~ 3 


0.0075 ± 0.0013 


A\/2Ar c /2, 


0.0107 


r/ -> p7 


U.Uol ± U.UUO 


n /in -L n no 
U.4U ± 


3B sin 0p 


U.4Z 


n' — > /ij'i/ 
r/ —> uj^ 


(ft i _i_ n o\ x iq-3 


n 11 + 01 


B[cos4>v sin0p±2r s sin0y cos^p] 


15 


p ± — > 7r ± 7 


0.068 ± 0.007 


0.22 ±0.01 


B 


0.21 


p° -► r?7 


0.036 ± 0.002 


0.45 ± 0.07 


3B cos 4>p 


0.46 


LJ — ► 7T°7 


0.71 ±0.04 


0.70 ±0.02 


3B cos 0y 


0.62 


(J — > 777 


(5.5 ± 1) x 10 


0.16 ± 0.01 


B [cos p cos <py — 2r s sm p sin 0y J 


n i i 

0.14 


— > 7T°7 


(5.8 ±0.5) x 10~ 3 


0.042 ± 0.002 


3.B sin 0v 


0.04 


-► 7?7 


0.056 ± 0.002 


0.21 ±0.01 


1? [2r s sin p cos 0y ± sin <p y cos p] 


0.20 


-> 7/7 


(5 ± 3) x 10~ 4 


0.30 ± 0.08 


B [sin v sin <pp — 2r s cos <^> y cos p] 


0.21 


- 7tV° 


< 0.23 ±0.01 


< 1.18 ±0.02 


3C sin <f)y 


0.69 


tp -> 7? c 7 


(1.1 ±0.4) x 10~ 3 


0.17 ±0.03 


4Br c 


0.19 


K*+ -► i^T+7 


0.050 ± 0.006 


0.252 ±0.015 


Br s [2J(r s )-J(l/r s )} 
(l+r s ) 


0.20 


K*° -► i^°7 


0.116 ±0.011 


0.39 ±0.02 


Br s [J(r s )+J(l/r s )] 

(l+r s ) 
Br c [-J(r c )+2J(l/r c )] 
1+r-c) 


0.34 


jD *+ _^ jD + 7 


< 0.0014 


< 0.14 ±0.10 


0.03 




< 0.8 

< 1.8 


< 3.4 

< 5.0 


2Br c [J(r c )+J(l/r c )] 

(1+Tc) 

Br s r c [J(r c /r s )+2J(r s /r c )] 
(r s +r c ) 


0.26 
0.04 



14 



TABLE II. Rare decays which are sensitive to the four mixing parameters 5 and e in (14) and 
(15). A, B and C are the same abbreviations as in Table I, and again r q = rh/m q stands for the 
constituent quark mass ratio. In writing down the ip amplitudes, we have set cos 4>v = 1 as a first 
approximation and ignored terms of order eS, which can be justified a posteriori. 



Process 


r ea; P(keV) 


\M exp \ (TeV^ 1 ) 




Mpw 


|M*^| (TeV^ 1 ) 


Vc -> P°P° 


(113 ±40) 


37 ±6 




—3C{5 C cos 4>p + 5' c sin^p) 


37 


r] c — > ujijj 


< 84 


< 32 




— 3C(<5 c cos 4>p + 5' c sm4>p) 


37 


Vc -> H 


94 ±30 


43 ±6 




3v / 2C(5 c sin0p — S' c cos <fip)r s 


43 


7/ c — > A A 


00 ± oU 


1 7 4- A 
1 ( ± 4 


3Cr s 
l+r s 


— (5 C cos 0p + S' c sin <f>p)J(r s ) 
+^(<5 c sin^p - 5^ cos <f> P )J(l/r s ) 


ID 


-0 — > 7T°7 


0.003 ± 0.001 


0.19 ±0.02 




-3Be c 


0.13 


— ► 


U.U 1 It U.UUo 


n Q9 4- n 


D 


-e c cos 0p - 2r s e' c sin (j) P + 4r c 5 c ] 


1 1 


ip — ► 7/7 


0.37 ±0.04 


2.3 ±0.1 


B 


— e c sin <fip + 2r s e' c cos 4>p + 4r c 5' c ] 


0.1 


V> -► 


0.37 ±0.02 


2.1 ±0.1 




-3Ce c 


2.1 


il> ^ rjuj 


0.14 ±0.02 


1.4 ±0.1 


—3C[e c cos <pp + V2r s e' c sin </>y sin <j)p] 


1.6 


tp — > r/0 


0.057 ± 0.008 


0.97 ±0.07 


—3C[e c cos 4>p siiKpy — y/2r s e' c sin 4>p] 


0.62 


^ — > r/u; 


0.015 ±0.002 


0.52 ±0.04 


—3C[e c sin <\>p — \^2r s e' c sin <j)y cos 4>p] 


1.4 


^ — > r/0 


0.029 ± 0.003 


0.80 ± 0.04 


— 3C[e c sin <pp sm<j)y + \f2r s e' c cos 4>p] 


0.68 


^ -► K°K*° 


0.37 ±0.03 


2.3 ±0.1 




6Cr s [e c J{r a )+e' c J(l/r s )] 
(l+rs) 


2.2 


A'+i?*- 


0.43 ± 0.07 


2.5 ±0.2 




6Cr a [e c J(r 3 )+e' c J(l/r s )} 


2.2 
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FIGURES 




. 5 



1 2 3 4 5 6 

r 

FIG. 1. Plot of the function J(r) from r = to 6. 
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